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The intermetallic FeSi exhibits an unusual tem- 
perature dependence in its electronic and mag- 
netic degrees of freedom, epitomized by the 
crossover from a low temperature non-magnetic 
semiconductor to a high temperature paramag- 
netic metal with a Curie- Weiss like susceptibility. 
Many proposals for this unconventional behav- 
ior have been advanced, yet a consensus remains 
elusive. Using realistic many-body calculations, 
we here reproduce the signatures of the metal- 
insulator crossover in various observables : the 
spectral function, the optical conductivity, the 
spin susceptibility, and the Seebeck coefficient. 
Validated by quantitative agreement with experi- 
ment, we then address the underlying microscopic 
picture. We propose a new scenario in which 
FeSi is a band-insulator at low temperatures and 
is metalized with increasing temperature through 
correlation induced incoherence. We explain that 
the emergent incoherence is linked to the unlock- 
ing of iron fluctuating moments which are almost 
temperature independent at short time scales. Fi- 
nally, we make explicit suggestions for improving 
the thermoelectric performance of FeSi based sys- 
tems. 

Iron based narrow gap semiconductors such as FeSi, 
FeSb2, or FeGa 3 show a pronounced resemblance to 
heavy fermion Kondo insulators in their charge and 
spin 0, [H degrees of freedom. Besides, these systems 
display a large thermopower 0, B BB El, heralding 
their potential use in solid state devices. There are two 
complementary approaches for explaining these unusual 
properties : On one hand it has been proposed that lat- 
tice degrees of freedom play a crucial role @, G3, 0] . On 
the other hand, electron-electron correlation effects have 
been invoked on the basis of both experimental results 




or a 



0, II] , as well as theoretical model studies 12 131 llj 
advocating in particular Hubbard physics 12t \ll 
spin fluctuations [IH , spin-state transitions la, |Tl 
thermally induced mixed valence [HI . 

Here, we go beyond modelistic approaches and investi- 
gate the effect of correlations on prototypical FeSi from 
the ab initio perspective. The key issues that we ad- 
dress are (1) can electronic Coulomb correlations alone 
quantitatively account for the signatures of the temper- 
ature induced crossover in various observables, and (2) 
what is the underlying microscopic origin of this behav- 
ior ? As a realistic many-body approach, we employ the 



combination of density functional theory and dynamical 
mean- field theory DFT+DMFT (for a review see e.g. Q) 
as implemented in Ref. @ . For the calculation of the See- 
beck coefficient we have extended our previous workfll| 
for DFT computations to include the DMFT self-energy, 
in a full orbital setup. For details see the supporting 
information. 

At low temperatures, iron silicide is a semiconductor 
with a gap A 50 - 60meV 0, 0, H3, with the resis- 
tivity [l|, [g and the magnetic susceptibility @ follow- 
ing activation laws. At 150-200K (i.e. at temperatures 
much smaller than A) a crossover to a (bad) metal is 
observed in transport [H, B , and optical spectroscopy 



HEEO-i!]. Moreover, FeSi displays a maximum in the 
susceptibility at 400K, followed by a Curie- Weiss like law 
Q. The energy scale over which spectral weight in the 
optical conductivity is transfered through the transition 
has been long debated [1, [l7 , 21 , 22 1 . Recent ellipsometry 
measurements 24[ showed that weight is moved over sev- 
eral eV - a common harbinger of correlation effects [25| . 
The Seebeck coefficient of FeSi peaks near 50K with a re- 
markable 700/iV/K, but is quickly suppressed at higher 
temperatures.when the unusual spin and charge proper- 
ties set in UEHzl. 
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FIG. 1: Local spectra and susceptibility. The theoretical 
local spectral function for different temperatures. uj m in traces 
the spectral minimum with respect to the Fermi level (lOx 
magnified), fnset : local spin susceptibility. 
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SPECTRAL PROPERTIES 

Fig. [T] shows our theoretical local spectral function of 
FeSi at various temperatures [orbital and momentum re- 
solved spectra can be found in the supporting informa- 
tion]. At low temperature, the spectrum is similar to that 
obtained within band-theory [16l4l8| . albeit with a gap 
renormalizcd by a factor of about 2, in agreement with 
photoemission spectroscopy (PES) (27, 28 1. [42| The elec- 
tronic excitation spectrum is thus band-like and coherent 
at low temperature. As was noted earlier 16-1^, the gap 
edges are very sharp in this regime, indicating a poten- 
tially large thermopower (29[. At higher temperatures, 
features broaden, and the system becomes a bad metal as 



found experimentally [27 



28| . We stress that this effect 
is largely beyond a mere temperature broadening of the 
electron distribution (Fermi) function 0, [27 , 28 1 . Due to 
the asymmetry of the spectrum, the position of the gap 
moves with temperature : In agreement with PES p^lici ] 
the minimum uj m i n of the spectrum [depicted in Fig. [I] 
starts out near the middle of the gap (as expected for a 
semiconductor [1 If ) . then first moves up with increasing 
temperature. Above 300K, uj m i n again approaches the 
Fermi level, as the asymmetry is reversed. 



OPTICAL SPECTROSCOPY 



compensation of excess carriers is achieved for the theo- 
retical curve. As is clear from Fig. [3[ AN(u>) does not 
vanish over the extended energy range plotted, hence a 
total compensation is not reached below the scale of the 
Coulomb repulsion of 5 eV. As a further assessment, we 
show, in Fig. Uta), the theoretical resistivity in compari- 
son with several experiments. 
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FIG. 2: Optical conductivity. Theoretical (top) and exper- 
imental [23 . |24( [ (bottom) optical conductivity as a function 
of frequency for various temperatures. 



Owing to a high precision, and the existence of sum- 
rule arguments, optical spectroscopy is a valuable tool 
for tracking the evolution of a system under change of 
external parameters. Transfers of spectral weight over 
scales related to the Coulomb interaction U rather than 
the gap A are usually considered a hallmark of correla- 
tion effects [25|] . Recent ellipsometry measurements [24| 
showed spectral weight transfers over several eV. In Fig. [2] 
we compare our realistic optical conductivities ct(uj, T) to 
experiment 23J, [24[ . Both, the overall magnitude and the 
crucial temperature evolution are well captured : Cooling 
down from 400K depletes spectral weight below 80meV, 
with only parts of it being transfered to energies just 
slightly above the gap. To analyze this in more detail, 
we follow Menzel et al. [24[ and plot in Fig. [3] the tem- 
perature difference AJV(w) = Nt x (w) — Nt 2 (w) of the ef- 
fective number of carriers Nt(oj) — 2 ™^¥ Jq du}'a(uj' , T) 
as a function of energy. An intersection with the x-axis 
corresponds to a full recovery of spectral weight as is im- 
posed by the f-sum rule. Our theoretical results quanti- 
tatively trace the experimental temperature dependence. 
We note that there are several isosbectic points in the 
optical conductivity, Fig. [2l which lead to extrema in 
AiV(w) in Fig. H The first peak in AJV is at 80meV, 
the scale of the semiconducting gap, above which spec- 
tral weight starts to pile up at low temperature in a(ui). 
The first minimum in AN(uj) occurs at the second isos- 
bectic point at around 0.1 8eV, up to where only ~ 35% 



Theory (580K - 77K) 
Experiment (200K - 7K) 




FIG. 3: Transfer of spectral weight. Difference of high 
and low temperature carrier density as obtained from the in- 
tegrated spectral weight [see text for details]. Experimental 
results from [24|, used temperatures as indicated. 



THERMOPOWER 

FeSi boasts a notably large thermopower at low tem- 
peratures [J H, 01 ■ Yet, it is compatible in magnitude 
with a Seebeck coefficient of purely electronic origin : For 
a band-like semiconductor in the regime ksT < A, the 
electron diffusive thermopower cannot exceed A/t(Ti|. 
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Prom this perspective, a larger gap favors the ther- 
mopower. While validated in FeSi, the above constraint 
is violated in the related compound FeSb 2 @, H, fTlj ] . 
Given the dominantly electronic picture, in conjunction 
with the band-like nature of FeSi at low temperature, the 
Seebeck coefficient in that regime can be accounted for 
by band theory : Indeed, it was shown that a slightly hole 
doped band-structure yields good agreement below 100K 
[171 ] . As shown in Fig. [5j we confirm this by supplement- 
ing band-theory with an effective mass of 2 and 0.001 
holes/Fe. While a dependence on stoichiometry is seen 
in experiments @, the tiny amounts of extra holes should 
be viewed as a means to alter the particle-hole asymme- 
try ll| rather than an effect of excess charge. As is well 
known, the particle-hole asymmetry plays a major role 
in determining the thcrmopower, since electron and hole 
contributions have opposite signs. The convention is that 
a negative Seebeck coefficient is dominated by electron 
transport, and a positive one by holes. It is interesting 
to note that for an insulator a large thermopower is ex- 
pected near thermoelectric particle/hole symmetry 1 lj . 
with a large sensitivity to the exact imbalance of carri- 
ers. Scanning through electron and hole doping, FeSi is 
experimentally indeed placed near such a boundary [30| . 

As expected, introducing an effective mass alone fails 
at higher temperatures. It has been conjectured (but 
not calculated) that this could be accounted for by ther- 
mal disorder via the electron- phonon coupling [ItI . [2o| . 
On the other hand, model studies indicated compatibility 
with the picture of Coulomb correlations (3l| . In Fig. [S 
we display the Seebeck coefficient obtained from our re- 
alistic many-body calculation. The agreement with ex- 
periments is very good. Notably, the theory captures 
the overall suppression of the thermopower above 100K, 
which comes from an enhanced conductivity due to the 
accumulation of incoherent weight at the Fermi level. In 
this sense incoherence is detrimental to semi-conductor 
based thermopower. Note however that correlation ef- 
fects can substantially enhance the Seebeck coefficient of 
correlated metals [see e.g. [32l|]. 

The temperature dependence of the Seebeck coefficient 
in Fig. [5] is connected to the moving of the chemical po- 
tential discussed above : As a function of rising temper- 
ature, the chemical potential first moves down, thereby 
increasing the electron contribution to the thermopower 



111 ]. At some temperature thermoelectric particle- hole 
symmetry is passed, and the Seebeck coefficient is domi- 
nated by electrons {S < 0), while at higher temperatures, 
the trend is reversed, and S becomes positive again. 

We note that the conversion efficiency of thermoelec- 
tric devices is measured by the so-called figure of merit 
ZT = S 2 oT/k. Here, the combination S 2 a of the See- 
beck coefficient 5" and the conductivity a is called the 
power factor (PF), and k is the thermal conductivity. 
The power factor of FeSi, displayed in Fig.@Jb), peaks at 
around 60K, where is reaches more than AO^iW / {K 2 cm), 



i.e. it reaches values comparable to state of the art Bi2Tc3 
(at its maximum value at 550K), and is for T > 60K 
larger than in FeSb2 that holds the overall record PF 
(realized at 12K)[3]. 



MICROSCOPIC INSIGHTS 

The quantitative agreement of our theoretical results 
with a large panoply of experimental data validates our 
approach, thus signaling the paramount influence of elec- 
tronic correlation effects. Drawing from the microscopic 
insights of our method, we now address the physical pic- 
ture underlying the intriguing properties of FeSi. This 
will in particular allow us to propose ways to improve 
the thermoclcctrical properties of FeSi. 



Crossover to the Metallic State 

Within our theoretical picture, the crossover to bad 
metallic behavior is not caused by a narrowing of the ex- 
citation gap [see the supporting information for details]. 
Instead, it is filled with incoherent weight that emerges 
with increasing temperature. Information about the co- 
herence of the one-particle excitations are encoded in the 
imaginary parts of the electron self-energy E. For the 
relevant orbital components we find 3£(u; = 0) ~ —aT 2 
with a = 1.9 ■ 10~ 4 mcV/K 2 . thus reaches the value 
of A/2 at around 400K, when only a pseudogap remains 
[see Fig. HJ. 

Complementarily, it was proposed that the arguably 
large electron-phonon coupling [7 , 2(J 24 1 causes the clo- 
sure of the gap via thermal induced atomic disorder 
17 . |20T |. In molecular dynamics simulations [2fj, the 
gap Ajjft was shown to vanish abruptly for tempera- 
tures of the order of T sa &. DFT /2, in contrast to the 
gradual transition that is observed in experiment and 
reproduced by our theory. We note that also in other 
systems with large electron-phonon coupling, spectral 
weight transfers arequantitatively accounted for by elec- 
tronic correlations 1331 . 



Strength of Correlations 

A recurring question in condensed matter physics 
is whether a system is well described in either an 
itinerant or a localized picture. Both the low effec- 
tive mass of FeSi (m*/m«1.5) and the rather high 
kinetic energy Eki n ~— 10.5eV of the iron states sig- 
nal a large degree of derealization. Indeed, other 
iron compounds show significantly higher masses and 
lower kinetic energies, e.g. the pnictides BaFe2As2 
{E kin [eV];(m*/m) xy )&{-7;3)&nd CaFe 2 As 2 (-8; 2.5) or 
the chalcogen FeTe (-6.5; 7) j2^|. Accordingly, FeSi is 
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an only moderately correlated, itinerant material. Yet, 
effects of correlation induced incoherence are essential 
for the crossover to the metallic state. This seeming 
contradiction is resolved by noting that all relevant en- 
ergy scales are of similar magnitude. Indeed : 0(A) re 
C(3S) re Q(T) re O(50meV), which leads low energy 
properties to be correlation dominated. 

The Spin State 



iron pnictides and chalcogenides[22|,[23j], and is quite dif- 
ferent from Hubbard physics in which the Hubbard U 
reduces charge fluctuations and the time-scale associated 
with them, while the time averaged local spin susceptibil- 
ity is strongly enhanced. In contrast, only the spin fluc- 
tuations but not the charge fluctuations are short-lived 
in FeSi, resulting in only a moderate enhancement of the 
(time averaged) local spin susceptibility, but a strongly 
enhanced short time (energy averaged) fluctuating mag- 
netic moment. 



A major signature of the unconventional behavior of 
FeSi is the non-monotonous uniform magnetic suscepti- 
bility @. The latter is closely mimicked by our local 
spin susceptibility xioc{u = 0) ~ ^ / cIt(S z (t)S z (0)) [see 
inset of Fig. Q], indicating compatibility with the pic- 
ture of electronic correlations. Interestingly there is a 
distinction of time scales : While xioc(u = 0) - the time- 
averaged response - displays a strong temperature de- 
pendence, the spin response at very short time scales, as 
probed by the observable lim r ^. (5'(' 7 ")5'(0)), is virtually 
temperature independent. From the latter, we obtain 
an effective moment M = \/S(S + l)g s re 3 (g s = 2). 
This is consistent with major contributions from effec- 
tive iron states with 5 = 1 [sec supporting information 
for details], and in agreement with M = 2.7 as obtained 
from fitting the experimental susceptibility Q to a Curie- 
Weiss law X = l f^-M 2 /{T~T c ) for T > 400K. We note 
that the local susceptibility [see inset of Fig. [T] is about 
one order of magnitude smaller than the experimental 
uniform susceptibility. Hence, the temperature induced 
fluctuating moment of the underlying spin state is, to a 
large extent, not local. This is corroborated by neutron 
experiments [IHl that find a significant magnetic scatter- 
ing at "ferromagnetic" reciprocal lattice vectors. 

Since the effective moment M is constant, FeSi thus 
does not undergo a spin-state transition as is found e.g. 
in MnO or LaCoC>3. There it originates from a competi- 
tion between the Hund's coupling J and the crystal field 
splittings. In FeSi, the splitting of low energy excitations 
is given by the gap A, which is much smaller than the 
Hund's coupling J, thereby favoring the high spin con- 
figuration. 

The preponderance of S=l states in particular implies 
that FeSi is not a singlet insulator as previously proposed 
IEEzI]. We further find that FeS i is in a mixed valence 
state. We obtain an iron valence of TV<j re 6.2, with a 
large variance SN = ((TV — (TV)) 2 ) re 0.93 [see support- 
ing information for details]. However, these numbers are 
insensitive to temperature, thus excluding a thermally 
induced mixed valence [3] ■ 

Interestingly, we find that our results are one order of 
magnitude more sensitive to the strength of the Hund's 
rule coupling J than to the Hubbard interaction U [for 
details see the supporting information]. This effect is 
reminiscent of the physics of Hund's metals as found in 



THE PHYSICAL PICTURE 

We are now in the position to elucidate the fundamen- 
tal picture of FeSi : (a) What is the physical origin of 
the crossover in the susceptibility ? (b) How is the latter 
linked to the emergence of incoherent spectral weight, i.e. 
what is the relation between the crossovers for the spin 
and charge degrees of freedom ? 

At low temperatures, FeSi is a conventional band-like 
semiconductor, i.e. it is approximately described by an 
effective one-particle Hamiltonian that can be diagonal- 
ized in momentum space. Excitations are well-defined 
(coherent). The system is in a high spin state, but spin 
excitations are gapped, and so the spin susceptibility is 
small as fluctuations at finite time scales are quenched. In 
this sense the spin degrees of freedom are inactive. With 
increasing temperature, however, a fluctuating moment 
develops at the iron sites. The emergence of such a lock- 
ing to the real-space lattice breaks down the momentum- 
space description : k is no longer a good quantum num- 
ber, hence excitations acquire a finite lifetime broaden- 
ing, and the spectrum becomes incoherent. 

Out theory goes beyond previous theoretical proposals, 
indeed it reconciles (and validates by including explicitly 
the Si degrees of freedom) two seminal model-based ap- 
proaches : The model of Fu and Doniach 12[ that ad- 
dressed the evolution of the one particle spectrum, and 
which can serve as a simplified cartoon of our realistic 
calculation, and the spin-fluctuation theory of Takahashi 
and Moriya 13| which is compatible with our results for 
the spin degrees of freedom. Here, we claim that the 
crossovers in the spin and charge response are actually 
intimately linked to each other. 

Due to the complicated multi-orbital nature of the 
problem, there arc however important differences from 
previous model calculations 12|, [l5|, [l8[ : In the realistic 



case the relevant control parameter is the Hund's cou- 
pling J rather than the Coulomb interaction U [see the 
supporting information]. Furthermore, we account for 
the particle/hole asymmetry. Indeed, for a symmetric 
model [H] , the Seebeck coefficient - a key feature of this 
system - is zero at all temperatures, and the chemical 
potential is pinned to its value at OK. 
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FIG. 4: Resistivity and Powerfactor. (a) Shown is a comparison of the theoretical resistivity of FeSi with several ex- 
serimental results[J, @, [I2> E3l- (b) the powerfactor of FeSi, theory in comparison with experimental data compiled from 
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FIG. 5: Thermopower. Theoretical Seebeck coefficient in 
comparison with experiment Q],E10|- Note the change in scale 
above 40//V/K. The line connecting the many-body theory 
results is a guide to the eye. The unconnected points at low 
temperature are obtained from an analysis of the asymptotic 
behavior of the many-body calculation. 



goal, we propose experimental studies of FeSi under ex- 
ternal pressure [39| , or by compressing the lattice using 
a suitable substrate. A particularly interesting option 
would be partial iso-electronic substitution of iron with 
ruthenium |4pj. Not only does the iso-structural RuSi 
have a larger gap than FeSi, and itself a notable See- 
beck coefficient above 100K[4lj, but also the fluctuat- 
ing moment physics that drives the metalization and the 
quenching of the thermopower in FeSi could be effec- 
tively reduced, as well as the thermal conductivity de- 
creased via the alloying. To our knowledge these avenues 
have not been pursued vis-a-vis their impact onto ther- 
moelectricity. More generally, the theory outlined here, 
describes the subtle interplay of electron-electron corre- 
lations and thermoelectricity, adding a new host of sys- 
tems amenable to theory assisted thermoelectric material 
design. 
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CONCLUSIONS & OUTLOOK 

To summarize, we have obtained a fundamental micro- 
scopic theory of iron silicide, the canonical example of a 
correlated insulator. This material features a competitive 
power factor in the temperature range of 50-100K. Our 
theory suggests concrete ways to improve that already 
remarkable thermoelectric performance. First, as is be- 
coming standard practice in thermoelectric development 
[jij , nano- or heterostructuring is needed to reduce ther- 
mal conductivity to convert the good power factor into 
a large figure of merit. More important yet, the power 
factor can be improved further by increasing the charge 
gap, through a reduction of the ratio of Hund's coupling 
over bandwidth, or an increase in the inter-atomic hy- 
bridizations that cause the band-gap. To achieve this 
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In wicn2k[l[ the Kohn-Sham orbitals are expressed in 
augmented plane waves (and local orbitals). Thus the 
above integral has contributions from within muffin tin 
spheres, interstitials and mixed terms. For details see 
Ref . [l5| . The spectral functions are given by 



METHOD : ELECTRONIC STRUCTURE 



For the density functional theory (DFT) part, we 
employ Wicn2k[l| within the generalized gradient 
approximation(GGA)0. The dynamical mean field the- 
ory (DMFT)0] impurity problem is solved using a hy- 
bridization expansion continuous time quantum montc 
carlo (ctqmc) method 0, Q. We use a projection based 
DFT+DMFT setup with full charge self-consistency, as 
implemented in Ref. 0] . For the Coulomb interaction and 
Hund's coupling, we use the value U = 5.0eV, J = 0.7eV, 
respectively, which have been found appropriate for iron 
based compounds 0]- For the dependence of results on J, 
see the discussion in the extended results section, below. 



METHOD : THERMOPOWER CALCULATIONS 

The electronic contribution to the thermopower along 
the a-direction, i.e. the diagonal element S a = S nn of 



the Seebeck tensor, can be written as (see e.g. Ref . [8j-|l2j ) 



Q " \e\ Ag 



(1) 



where -Aq/i arc ^ nc coefficients that relate (within linear 
response theory) the perturbed charge/heat current to 
the external perturbing field. They can be expressed as 



A" 



dui 



HaH (2) 



with the Fermi level [i and the Fermi function / M . When 
neglecting vertex corrections in the DMFT spirit [l3, 14 1, 
the transport kernel can be written as 



HaCw) = £Tr[u Q (k)A(k,w)i; Q (k),4(k,w)] (3) 



where the transition matrix elements (Fermi velocities) 
i> Q (k) and the spectral functions A(k,uj) are matrices in 
orbital space. We choose to evaluate the trace in 3 in 
the Kohn-Sham basis {^'kri(r)}. Accordingly, the Fermi 
velocities - that in the dipole approximation are matrix 
elements of the momentum operator V - are computed 
as 

1 



v™ (k) = -(kn\V a \kn' 



= / d 3 r*^ n (r)V«* kll ,(r) 

m J 



A(k,u) = -— [G(k,«)-G+(k,«)] (5) 
with the Greens function G in the Kohn-Sham basis 



G~ 



(k,w) = (lu + /i - e km )6r, 



,(k,w) 
(6) 



where ekm are the Kohn-Sham energies. S = H — Ejjc 
is the double counting corrected self-energy obtained by 
embedding the local impurity self-energy S T (w) of atom 
t into the solid by (see Ref. Q for further details) 

£n»'(k,w) = Y, PZ(nn',LL')Xl L ,(uj) (7) 

tLL' 

where L, L' index the impurity space, in our case referring 
to spherical harmonics L = (l,m). Here, P is the projec- 
tion operator that defines the Hilbcrt space of the corre- 
lated orbitals, as well as the embedding procedure @. In 
our DFT+DMFT calculation, we employ full charge self- 
consistency, whereby not only the spectral functions are 
influenced by the correlation effects (via the self-energy), 
but also the transition matrix elements are altered by 
means of a self-consistent charge density. 



FURTHER DETAILS OF THE RESULTS 

Band- structure. FeSi crystallizes in the B20 structure 
with the cubic space group P2i3, and four iron atoms 
per unit cell. The ligand arrangement splits the 3d or- 
bitals into 1/2/2 for which we find the corresponding 
major characters z 2 /x 2 -y 2 + xy/xz + yz, with the latter 
two multiplets separated by an avoided crossing. With 
a nominal iron valence of 6, FeSi is an insulator within 
band theory (l6l - [l8j . Owing to the hybridization nature 
of the gap, and the absence of bandwidth narrowing ef- 
fects -and contrary to most other semi-conductors-, DFT 
calculations overestimate the gap by a factor of about 2 : 
A DFT « O.lleVlIHii. 



(4) 



Orbital and momentum resolved Spectra and Influ- 
ence of Incoherence. In Fig. [5] (a) we display the 
DFT+DMFT spectral function, resolved into orbital 
characters of the Fe 3d manifold (different line types) 
and for two temperatures (red and blue). The gap (or 
pseudo gap, depending on the temperature) is formed 
mainly by orbitals of x 2 — y 2 and xy character on the 
valence, and xz and yz on the conduction side. Besides 
the moving of the chemical potential [see discussion in 
the main text], the spectral weight at the Fermi level 
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is increased at larger temperatures by effects of incoher- 
ence rather than a closing of the gap though a shifting of 
excitations. In the following, we shall justify this claim. 

The impurity Greens function (of atom r and orbital 
L = (I, m)) is given by 

GI l (uj) = [c-K^ + AElH + S^H)]- 1 (8) 

where E imp are the impurity levels, A(w) the hybridiza- 
tion function, and £(w) the self-energy, and we have as- 
sumed all quantities to be diagonal in the local basis (for 
details see e.g. Q). A good idea about the major fea- 
tures in the excitation spectrum can be gained (see e.g. 
Ref.[lj|) by analyzing the poles of the Greens function 
in the absence of lifetime effects, i.e. by 



L 



(9) 



which is shown in Fig. [Fj] (b) where the intersections with 
the black line mark the graphical solutions for the above 
equation. As expected, poles in the occupied part are of 
x 2 — y 2 and xy character, whereas as for uj > solution 
can be found for the xz and yz components. 

The interesting observation here is that those poles 
slightly move further apart when the temperature is 
raised, i.e. the gap in the quasi-particle band-structure 
actually becomes larger. Yet at the same time, as is 
evident from Fig. [6] (c), the imaginary part of the self- 
energy increases drastically. The values of SjE are seem- 
ingly small, but they are comparable in magnitude to the 
experimental charge gap (50 — 55meV). Fig. [6jd) shows 
that the values at the Fermi level of the relevant orbital 
components of follow up to 400K a T 2 law with a co- 
efficient a = —1.9 • 10~ 4 meV/K 2 . Assuming Lorentzian 
line shapes, the combined half-width-at-half-maximum of 
the conduction and valence bands reaches the size of the 
gap below 400K : 23£(u; = 0,T = 363K) = 50meV. 
Yet, we note that the filling of the gap is gradual and al- 
ready at lower temperatures, there appears finite incoher- 
ent weight at the Fermi level. This is different from the 
closure of the gap in the picture of thermal disorder poj. 
mentioned in the manuscript, where the bridging of the 
gap occurs rather abruptly for T as A/2. Also, due to 
the particle/hole asymmetry, the chemical potential falls 
into the middle of the gap only at T=0 [see also the main 
text, and Fig.l there]. This further lowers the temper- 
ature at which sizable incoherent weight appears at low 
energy. 

For accessing lower temperatures in the Seebeck coef- 
ficient than is possible for the Monte Carlo to reach (see 
the main manuscript, and, there, the unconnected points 
in Fig. 4), we extrapolate the imaginary parts of the self- 
energy of the run at the lowest temperature following the 
T 2 behavior described above. 

In Fig. [71 we further show the total spectral function 
momentum-resolved along various symmetry lines, and 
for two temperatures. As was evident already from the 



local spectra, excitations broaden significantly with tem- 
perature. Also seen is the slight moving apart of the peak 
intensities at higher T, i.e. the widening of the gap as 
given by the poles of the Greens function discussed above. 
While the incoherent spectral weight at the Fermi level is 
not discernible on the used color scale, it is evident that 
the crossover to a (bad) metallic state is not caused by 
the introduction of a "band" crossing. 

Charge Carrier Concentration As a reference, we fur- 
ther provide, in Fig. [5J the temperature dependence of 
the carrier concentration, extracted from our data by two 
different means : 

(1) from the spectral function. Using n = 
V So° dujf(uj)A(uj) and p = ± J ^ duf(-u)A(u) with 
the unit-cell volume V and the local spectral function 
A((jj) we compute the number of electrons and holes, re- 
spectively. We find that the net carriers are holes at all 
temperatures (p — n > 0). This is to be contrasted with 
the Seebeck coefficient that from low to high tempera- 
tures changes from p-type to n-type and back to p-type, 
demonstrating that particle-hole symmetry for different 
response functions can be distinct. 

(2) from the restricted sum-rule of the optical conductiv- 
ity. Using J Q diua(uj) = §^n e ff with a gap A = GOmeV 
and the orbital averaged effective mass m* = 1.5m, we 
find values for the number of effective carriers that are 
in good agreement with those obtained from the one- 
particle spectrum. 

We however note that, in particular at low tempera- 
tures, extrinsic effects that are not included in our calcu- 
lation can potentially have a large effect on the number 
of carriers. Indeed the sign of the Hall coefficient seems 
to be notably sample dependent. 

Dependence of Results on the Hund's coupling. In 
FeSi there is no direct competition between the Hund's 
coupling J and the crystal field splittings, as e.g. manifest 
from the energy difference between the poles of the im- 
purity Greens function (see above). The latter is about 
0.3eV, while the Hund's coupling used for our calcula- 
tions is more than twice as large : J = 0.7cV. Still, there 
is a notable dependence on the value of J (see also the 
recent Ref. 2l|) : In Fig. [5] (a) are displaced the effective 



masses to* /to as obtained from the slope of the local self- 
energy at the Fermi level, as well as (b) the imaginary 
part of the self-energy at the Fermi level. 

As was found previously for other iron based 
as well as 4d systems [24|. the effec- 

The charge 
The 



materials [22L 125 

tive masses increase significantly with J. 
gap as shown in Fig. [§] (c) shrinks accordingly, 
general mechanism and system atics b ehind this are in- 
vestigated e.g. in Ref. [25 



24, 24-271, and can be ex- 



plained by resorting to a description in terms of a self- 
consistent Kondo problem [2 5|. A large Hund's rule cou- 
pling constraints electrons to same spin, different orbital 
states, thereby "orbitally blocking" [23j the Kondo inter- 
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action, which in turn decreases the coherence tempera- 
ture exponentiallvpBj . As a consequence, also the imagi- 
nary part of the self-energy, as shown in Fig.|H](b), grows 
significantly in magnitude when increasing J. For our 
main calculations, we use the value J = 0.7eV which is 
found to be typical for iron based compounds 0. 

Spin state. The figures Fig. flOta) and (b) contain 
the histograms of the ctqmc[U for two different tempera- 
tures : Shown is the probability distribution of the many- 
body wave-function onto the eigenstates of the effective 
iron atom representing the 3d orbitals, decomposed into 
the number of particles N and the spin state S. The main 
conclusions that we can draw from these results are 

• the decomposition onto N and S is almost inde- 
pendent of temperature, but strong fluctuations at 
short time-scales are present. 

• the system is strongly mixed valence, since 
the charge variance is of order 1 : SN = 
((N- (N)) 2 ) w0.93. 

• the spin state S=l has the overall largest proba- 
bility at all temperatures, yielding an effective mo- 
ment M = y/S(S + Vjg s w 3 (g s = 2), with vari- 
ance SS = ((S - (S)) 2 ) « 0.33. 

Thus, while the variances of S and N are large, they 
do not evolve with temperature. This is in strong con- 
trast to systems like MnO, or LaCoC>3 in which spin-state 
transitions occur 28l- 3l| . 

We note that without averaging over states with 
the same S, the single most probable state has 
(iV,S)=(7,1.5), closely followed by (N,S)=(Q,2). This 
is similar to other iron systems like e.g. the pnictidc 
BaFe9As9[23|. However, in BaFe2As2 the most promi- 
nent states have a probability that is larger by a factor 
of 4 than most other states [231 s while the distribution of 
states of FeSi is remarkably uniform, as indicated by the 
character of the state averaged S. 
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FIG. 6: (a)-(c) Spectral function, real, and imaginary parts 
of the DMFT self-energy resolved into orbital components are 
shown for the two temperatures T=116K (blue) and T=387K 
(red), (d) imaginary parts of the self-energy at zero frequency 
as a function of temperature. 
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FIG. 8: Temperature dependence of the carrier con- 
centration. Density of carriers extracted from the one- 
particle spectrum (red) and from the optical conductivity 
(green). See text for details. 
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FIG. 9: Dependence of the effective masses (left), the imaginary parts of the self-energy (middle) and the spectrum (right) on 
the Hund's rule J, for U = 5.0eV and T = 116K. 
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FIG. 10: Probability of atomic states of the DMFT impurity, decomposed in number of particles N (x-axis) and spin-state S 
(colors). Results for (a) T = 116K and (b) T = 1161AT. 



